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Quadratic Riccati Differential 
Equation 
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Abstract 


In this paper, we introduce an efficient method for solving the quadratic 
Riccati differential equation. In this technique, combination of Laplace trans- 
form and new homotopy perturbation methods (LTNHPM) are considered 
as an algorithm to the exact solution of the nonlinear Riccati equation. Un- 
like the previous approach for this problem, so-called NHPM, the present 
method, does not need the initial approximation to be defined as a power 
series. Four examples in different cases are given to demonstrate simplicity 
and efficiency of the proposed method. 


Keywords: Riccati differential equation; Laplace transform method; NHPM; 
LTNHPM. 


1 Introduction 


The quadraic Riccati differential equation, deriving its name from Jacopo 
Francesco, Count Riccati (1676-1754). These kinds of differential equations 
are a class of nonlinear differential equations of much importance, and play 
a significant role in many fields of applied science [1]. At an early stage, the 
occurrence of such differential equations in the study of Bessel functions led 
to its appearance in many related applications and, to the present time, the 
literature on the Riccati equation has been extensive. For several reasons, a 
Riccati equation comprises of a highly significant class of nonlinear ordinary 
differential equations. Firstly, this equation is closely related to ordinary 
linear homogeneous differential equation of the second order. Secondly, the 
solution of Riccati equation possesses a very particular structure in that the 
general solution is a fractional linear function of the constant of integra- 
tion. Thirdly, the solution of an Ricaati equation is involved in the reduction 
of nth-order linear homogeneous ordinary differential equations. Fourthly, a 
one-dimensional Schrodinger equation is closely related to Riccati differential 
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equation. Solitary wave solutions of a nonlinear partial differential equation 
can be expressed as a polynomial in two elementary functions satisfying a pro- 
jective Riccati equation [3]. Moreover, such types of problems also arise in the 
optimal control literature, and can be derived in solving a second-order lin- 
ear ordinary differential equation with constant coefficients [2]. In conformity 
with the general study of differential equations, much of the early works were 
concerned with the study of particular classes of Riccati differential equation 
with the aim of determining the solution in finite form. Many mathemati- 
cians such as Jame Bernouli, Jahn Bernouli (1667-1748), Leonhard Euler 
(1707-1783), Jean-le-Rondd’Alembert (1717-1783), and Adrian Marine Leg- 
endre (1752-1833) contributed to the study of such differential equations [1]. 
Deriving analytical solution for Riccati equation in an explicit form seems to 
be unlikely except for certain special situations. Of course, having known its 
one particular solution, its general solution can be easily derived. Therefore, 
one has to go for numerical techniques or approximate approaches for getting 
its solution. Recently, Adomian’s decomposition method has been proposed 
for solving Riccati differential equations [7, 8]. HPM was introduced by He 
4| and has been already used by many mathematicians and engineers to 
solve various functional equations [9, 10, 5, 6, 12, 11]. Abbasbandy solved a 
Riccati differential equation using He’s VIM, homotopy perturbation method 
HPM) and iterated He’s HPM and compared the accuracy of the obtained 
solution to that derived by Adomian decomposition method [13, 14, 16]. 
Moreover, Homotopy analysis method (HAM) and a piecewise variational it- 
eration method (VIM) are proposed for solving Riccati differential equations 
16]. In [15], Liao has shown that HPM equations are equivalent to HAM 
equations when i = —1. Aminikhah and Hemmatnezhad [6] proposed a new 
homotopy perturbation method (NHPM) to obtain the approximate solution 
of ordinary Riccati differential equation. In this work, we present the solu- 
tion of Riccati equation by combination of placeLaplace transform and new 
homotopy perturbation methods. An important property of the proposed 
method, which is clearly demonstrated in examples, is that spectral accuracy 
is accessible in solving specific nonlinear Riccati differential equations which 
have analytic solution functions. 


2 LTNHPM for quadratic Ricatti equation 


Consider the nonlinear Riccati differential equation as the following form 


fw) = A(t) + B(t)u(t) + C(Hu2(), O< t <T, 4) 


(0) = 


where A(t), B(t) and C(t) are continuous and qa is an arbitrary constant. By 
the new homotopy technique, we construct a homotopy U : 2 x [0,1] —-R, 
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which satisfies 
H(U(t),p) = U'(t) —uo(t) +puo(t) —p[A(t) + BU (t) + C(£)U?(t)] = 0, (2) 


where p € [0,1] is an embedding parameter, uo(t) is an initial approximation 
of solution of equation (1). Clearly, we have from equation (2) 


H(U(t),0) = U"(t) — uo(t) = 0, (3) 
H(U(t),1) = U'(t) — A) — B)U(t) — C()U?(d) = 0 (4) 
By applying Laplace transform on both sides of (2), we have 


L {U'(t) — uo(t) + puo(t) — p[A(t) + B(t)U(t) + C(HU7(t)]} =0 (5) 


Using the differential property of Laplace transform we have 
sL{U(t)}—U (0) = L {uo(t) — puo(t) + p[A@) + BHU) + CHU? (H)]} (6) 


or 


L{U(t)} = : {U(0) + L {uo(t) — puo(t) + p[A(t) + B)U(t) + CHU? (t))} } 
(7) 


By applying inverse Laplace transform on both sides of (7), we have 


Uli) = 1-40 (0) + Lwolt) puo(t) + p[A(t) + BU (t) + C(HU*(t)]}}} 
(8) 


According to the HPM, we use the embedding parameter p as a small pa- 
rameter, and assume that the solutions of equation (8) can be represented as 
a power series in p as U(t) = 0°) p"Un . Now let us write the Eq. (8) in 
the following form 


> P"Un(t) {Yoo +L{uo(t) — puo(t) 
n=0 


+ p[At) + BC) D> p"Ualt) + CD »"Ualt))?] }) (9) 


Comparing coefficients of terms with identical powers of p leads to 


ps Uld= LAs ( (U(0) +L {uo(t)})} 

pi: Ur(x) =L-* {1 (L {up(t) — A) — BDU Os } 

pls Uj(a) =L-1 {2 (L {BU ;-1(8) + CO) HAZ) Ue (ja “0})} 
re 


(10) 
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Suppose that the initial approximation has the form U(0) = uo(t) =a, 
therefore the exact solution may be obtained as following 


u(t) = lim U(¢) = Uo(t) + Ui(t) + -- (11) 


pol 


3 Examples 


Example 1. Consider the following quadratic Riccati differential equation 
taken from [1] 


u'(t) = 16t? — 5 + 8tu(t) + u?(t), (12) 
u(0) = 1 
The exact solution of above equation was found to be of the form 
u(t) = 1 — 4¢. (13) 


To solve equation (12) by the LTNHPM, we construct the following homotopy 
U'(t) — ug(t) + p [uo(t) + 5 — 16¢7 — 8tU(t) — U7(t)] =0 (14) 
Applying Laplace transform on both sides of (14), we have 
L {U'(t) — uo(t) + p [uo(t) +5 — 16t? — 8tU(t) — U?(t)] } =0 (15) 
Using the differential property of Laplace transform we have 
sL{U(t)} — U(0) = L {uo(t) — p [uo(t) + 5 — 16¢” — 8tU(t) — U7(t)]} (16) 
or 


1 
L{U(t)} = . {U(0) + L {uo(t) — p [uo(t) + 5 — 16¢? — 8tU(t) — U*(t)] }} 
(17) 
By applying inverse Laplace transform on both sides of (17), we have 


1 
CC ee { (U(0) + L {uo(t) — p [uo(t) + 5 — 1647 — 8tU(t) — U7(t)]|}) 
8 
(18) 
Suppose the solution of equation (18) to have the following form 
U(t) = Uo(t) + pUr(t) + p°Ua(t) +++, (19) 


where U;(t) are unknown functions which should be determined. Substituting 
equation (19) into equation (18), collecting the same powers of p and equating 
each coefficient of p to zero, results in 
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po: Up(t) =L* {5 (UO) + L {uo(t)})} 

(x) =L~! {—4 (L {uo(t) +5 — 16? — 8tUo(t) — UB(t)})} 
2 (L {-80U;-1(¢) + se Uj, (t)Uj—K-1 

(20) 

Assuming ug(t) = U(0) = 1, and solving the above equation for U;(t), 7 = 

0,1,--- leads to the result 


Uo(t) =1+¢, 
Ui (t) = (5? + 3t—3), 
Uo(t) = & (10¢3 + 102? — 8t — 3), 
U3(t) = © (425¢4 + 59523 — 3991? — 399¢ + 63) , 
(t) = 4 (38750¢° + 6975024 — 36000¢3 — 72135¢? + 7938¢ + 8505) , 


Therefore we gain the solution of Eq. (12) as 
u(t) = Uo(t) + Ur(t) + Us(t) +--+» = 1-4 
which is exact solution. 


Example 2. Consider the following quadratic Riccati differential equation 
taken from [13, 7, 14, 6, 8] 


pci 


The exact solution of above equation was found to be of the form 


1 Jf2-1 
u(t) =1+ V2tanh | V2t + = lo 22 
(t) 5 log ( aa ‘) (22) 
The Taylor expansion of u(t) about t = 0 gives 
1 t t 53 71 
=t+t?+4—— x1 a 4 t+ a Ph +-.-. (2 
uo) gf 3g’ ane a5 © B55 29) 


To solve equation (21), by the LTNHPM, we construct the following homo- 


topy 
U'(t) = uo(t) — p [uo(t) — 1 — 2U(t) + U?(0)] (24) 


Applying Laplace transform, we have 
L {U'(t) — uo(t) + puo(t) — p[l + 2U(t) — U7 (t)]} =0 (25) 


Using the differential property of Laplace transform we have 
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sL{U(t)} — U(0) = L {uo(t) — p[uo(t)-1-2U(t)+U7H]} (26) 
or 


1 
L{U(t)} = = {U(0) +L {uo(t) — p [wo(t) — 1 — 2U(t )+U°H)]}} (27) 
By applying inverse Laplace transform on both sides of (27), we have 
1 
U() <1 {£ (U(O) +L {uolt) ~p [ualt) - 1-20 + 0(0]})} (28) 
Suppose the solution of equation (28) to have the following form 


U(t) = Uo(t) + pUi(t) + p?Ua(t) +---, (29) 


where U;(t) are unknown functions which should be determined. Substitut- 
ing equation (29) into equation (28), collecting the same powers of p , and 
equating each coefficient of p to zero, results in 


Cae (U(0) pane 
one 11 (L ea 
3 r U;(2)= Le 142 (1 (20; 


ith 2 U;(t) U;-1-1@})}, 9=2)3,0. 

(30) 
Assuming u(t) = U(0) = 0 , and solving the above equation for U;(t), 7 = 
0,1,--- leads to the result 


Say 


Uo(t) = 0, 

Ui (t) _ t, 

Ua(t) = t?, 

U3(t) = a (31) 
Ua(t) =— 

Us(t) = SAGs 


Therefore we gain the solution of equation (21) as 


ey Ol) ce ae ae 7 46 53 47 (32) 
=t+P+48—-h4- Ce - 24+ Be +.-., 


and this in the limit of infinitely many terms, yields the exact solution of (21). 


Example 3. Consider the following quadratic Riccati differential equation 
taken from [1] 
u(t) = . — e%! + 2e*tu(t) — efu?(t), 


0)= (33) 
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The exact solution of this equation is u(t) =e’ . We construct the following 
homotopy 


U'(t) — uo(t) + p | uo(t) — e?* + e8* — 2eP'U(t) + S> eP'U7(t)| =0 (34) 
n=0 
Applying Laplace transform on both sides of (34), we have 
L {or — uo(t) +p | uo(t) — e?? + e8P* — 2e??*U(t) + S- e'U?(t) =0 
n=0 
(35) 
Using the differential property of Laplace transform we have 
sL{U(t)} — U(0) =f sa - p[uo(t) — ert 
462?! — Qe*P U7 (t) + S- er"'U(1)| (36) 
n=0 


or 
L{U(t)} Huw a3 L{wo(t) 


—p [uo (t) — oP + e8Pt — 26?tU(t) + > er"'U(1)| }} (37) 
n=0 


By applying inverse Laplace transform on both sides of (37) and using the 


Taylor series of e* = >> 9 fod)" , we have 
u(t) U4 (uy) + Lf u(t) - p[uo(®) - Soe" 
8 ae n! 
+p OP aS py 4 oh ve]})b os) 
err n! =a n! ae n! 


Suppose the solution of U can be expanded into infinite series as U(t) = 
ye P"Un(t) . Substituting U(t) into equation (38), collecting the same 
powers of p , and equating each coefficient of p to zero, results in 
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pr Gy aL Aa aa )})} 

pi: U(x) =L71{-< : uo(t (1) — 2Uo + US})} 

p?: Ux(e) =L“! L {2¢ ie QU, + tUs + 2UoU1})} 

p?: U3(z) ne Se ae AtU, — 2U2 — 4t2Up + 2tUoUr (39) 


it BU U2 + 4t?UG})} 


Assuming uo(t) = U(0) = 1, and solving the above equation for U;(t), 7 = 
0,1,--- leads to the result 


Uo(t) =1 +t, 

Ui (t) = “> 

U2(t) = 13,(- 6 —8t+t?), (40) 
U3(t) = — sty (840 — 4620t + 504¢? — 840t9 — 240¢4 + 35°) 


Therefore we gain the solution of equation (33) as 
u(t) _ Uo(t) t U,(E) t U3(E) fenaee 
Slttt gt get: (41) 
=n a 


which is exact solution. 
Example 4. Consider the following quadratic Riccati differential equation 
taken from [13, 7, 14, 6, 8] 


u'(t) = —u(t) + u(t), 

{0) =} ne 

The exact solution of above equation was found to be of the form 

-t 
e€ 
t) = 43 
u(t) = (43) 
The Taylor expansion of u(t) about t = 0 gives 
1 1 1 1 17 31 691 

i= t ie a i v4 pra... 
ue) 2 A 48 480 Le 80640. 1451520 319334400 (44) 


To solve equation (43) by the LTNHPM, we construct the following homo- 
topy: 
U'(t) = uo(t) — p [uo(t) + Ut) — U*(0)] (45) 


By applying Laplace transform we get 


L{U(t)} = : {U(0) +L {uo(t) — p [uo(t) + U(t) — U*(t)]}} (48) 
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Using inverse Laplace transform on both sides of (46), we have 


1 


UG) =L>4 {; (U(0) + L {uo(t) — p [uo(t) + U(t) — U7 (t)] )} (47) 


Suppose the solution of equation (47) to have the following form 
U(t) = Uo(t) + pU i (t) + p?Ua(t) + ---, (48) 


Substituting equation (48) into equation (47), collecting the same powers of 
p, and equating each coefficient of p to zero, results in 


p?: Uo(t) = aa (U(0) ae we 
pi: vile) = {2 (L {up ce o(t)})} 
p> U;(x)=L7 1/2 (L{-u, Uz, (t)U;—K- i) })}, i oO ee 


(49) 
Assuming uo(t) = U(0) = 4 , and solving the above equation for U;(t), 7 = 
0,1,--- leads to the result 


Uo(t) = 3. +¢), 

Ui (t) = 7p (- 9+ #7), 

U2(t) = 00, (-15 +27), 

U3(t) = 504 1 (945 — 387¢? + 17¢*) , (50) 
Uas(t) = Tap 5 (5103 — 918t? + 31¢4) , 

Us(t) = qoegaqp (— 280605 + 2272052? — 2557514 + 691°) , 


Therefore we gain the solution of equation (42) as 


u(t) = Up(t) + Ur(t) + Us(t) + --- 
ee ees ea ae WO i Ne ere 31494 __ 691 
2 4° 48 780 80640 1451520 319334400 


fil 


(51) 
this in the limit of infinitely many terms, yields the exact solution of (42). 


4 Conclusion 


In the present work, we proposed a combination of Laplace transform method 
and homotopy perturbation method to solve nonlinear Riccati differential 
equation. The new method developed in the current paper was tested on 
several examples. The obtained results show that these approaches can solve 
the problem effectively. Unlike the previous approach implemented by the 
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present authors, so-called NHPM, the present technique, does not need the 
initial approximation to be defined as a power series. In the NHPM we reach 
to a set of recurrent differential equations which must be solved consecutively 
to give the approximate solution of the problem. Sometimes we have to do 
many computations in order to reach to the higher orders of approximation 
with acceptable accuracy. Also as an advantage of the LTNHPM over decom- 
position procedure of Adomian, the former method provides the solution of 
the problem without calculating Adomian’s polynomials. The advantage of 
the LTNHPM over numerical methods (finite difference, finite element, ...) 
and ADM is that it solves the problem without any need to discretization 
of the variables. The Computations finally lead to a set of nonlinear equa- 
tions with one unspecified value in each equation. The results show that the 
LTNHPM is an effective mathematical tool which can play a very important 
role in nonlinear sciences. 
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